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Statistical Characterization and Mitigation of NLOS 
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Abstract — In this paper the problem of the joint statis- 
tical characterization of the NLOS bias and of the most 
representative features of LOS/NLOS UWB waveforms 
is investigated. In addition, the performance of various 
maximum-likelihood (ML) estimators for joint localization 
and NLOS bias mitigation is assessed. Our numerical 
results evidence that: a) the accuracy of all the considered 
estimators is appreciably affected by the LOS/NLOS 
conditions of the propagation environment; b) a statistical 
knowledge of multiple signal features can be exploited to 
mitigate the NLOS bias, so reducing the overall localization 
error. 

Index Terms — Ultra Wide Band (UWB), Radiolocaliza- 
tion. Bias mitigation. 



I. Introduction 

Wireless localization in harsh communication 
environments (e.g., in a building where wireless 
nodes are separated by concrete walls and other 
obstacles) can be appreciably affected by direct path 
attenuation and NLOS conditions. In principle, the 
effects of NLOS errors can be mitigated adopting 
techniques for detecting LOS/NLOS conditions or 
algorithms for estimating the errors themselves (i.e., 
the bias due to obstacles). These approaches can 
potentially improve the overall accuracy and are ex- 
pected to satisfy certain requirements; in particular, 
they should be robust against changes in the prop- 
agation environment, stable (i.e., they should never 
amplify NLOS errors) and should be able to exploit 
all the available data (e.g., received waveforms and 
a priori information). 

Various solutions for NLOS error mitigation in 
UWB environments are available in the technical 
literature El, [|2||, [|3|, H. In particular, a simple 
deterministic model, dubbed wall extra delay, is 
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proposed in [[T]| to estimate the bias introduced 
by walls. A non-parametric support vector ma- 
chine is employed in [2] for joint bias mitigation 
and channel status detection; this approach exploits 
multiple features extracted from received signals 
in a non-statistical fashion. A few classification 
algorithms for LOS/NLOS detection are compared 
in [3 1, where it is shown that the best solution 
is offered by a statistical strategy based on the 
joint probability density function (pdf) of the delay 
spread and the kurtosis extracted from the received 
signals. Finally, in [4] statistical models for the 
time of arrival (TOA), the received signal strength 
(RSS) and the root mean square delay spread (RDS) 
are developed and an iterative estimator for bias 
mitigation is devised. 

The contribution of this paper is twofold. In fact, 
first of all, the problem of joint statistical modeling 
of multiple features extracted from a database of 
waveforms acquired in a TOA-based localization 
system is investigated. Note that, as far as we know, 
in the technical literature only univariate models for 
bias mitigation have been proposed until now (e.g., 
see flU, [|4|| and dSj]). The use of multiple signal 
features in UWB localization systems has been 
investigated in [[S] for channel state detection only 
and in Q, where, however, the considered features 
(namely, the kurtosis, the mean excess delay and the 
delay spread) have been modelled as independent 
random variables. The second contribution offered 
by this manuscript is represented by a performance 
comparison of various maximum-likelihood (ML) 
estimators for TOA-based localization. In particular, 
unlike other papers (e.g., see |2l, ||3l) we illustrate 
some numerical results referring to the accuracy 
of different localization strategies, rather than to 
capability of the bias removal on a single radio link. 

The remaining part of this paper is organized as 
follows. In Section |ll] some information about our 
UWB experimental campaign and about the features 
extracted from the acquired data are provided. In 
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Section III some estimation algorithms for UWB 
radiolocalization are described, whereas their per- 
formance is compared in Section [TV] Finally, some 
conclusions are given in Section |V[ 

II. Experimental setup 
A. Measurement arrangement 

A measurement campaign was conducted by our 
research group in the second half of 2010 in order 
to assess the impact of NLOS bias errors on a 
real world UWB system for localization; all the 
measured data were acquired by means of two 
FCC-compliant PulsON220 radios commercialised 
by TimeDomain and were collected in a database. 
Such devices are equipped with omnidirectional 
antennas, are characterized by a -10 dB bandwidth 
and a central frequency equal to 3.2 GHz and 
4.7 GHz, respectively, and perform two-way TOA 
ranging; they also allow to store digitised received 
waveforms (a sampling frequency of 24.2 GHz and 
14 bits per sample are used). 

It is worth mentioning that various databases 
providing a collection of sampled UWB waveforms 
acquired in experimental campaigns and useful for 
assessing the performance of localization algorithms 
are already available (e.g., see 10, [HI). However, 
our database has been specifically generated to 
assess the correlation between the NLOS bias error 
and various features extracted from the received 
signals, as it will become clearer in the next Para- 
graph. Our measurement campaign consisted of two 
phases. First, the transmitter was placed in a given 
room (room A in Fig. [T|) and the receiver in an 
adjacent room (room B in Fig. [T]) separated from 
room A by a wall having thickness tyjaii = 32 cm 
(NLOS condition); in addition, the transmit antenna 
was kept fixed, whereas the receive antenna was 
placed in A^^^g^ = 174 distinct vertices of a dense 
square grid (the distance between a couple of nearest 
vertices was equal to 21 cm). In the second phase 
of our measurement campaign both the transmitter 
and the receiver were placed in room B; then, the 
transmit antenna was kept fixed, whereas the receive 
antenna was moved on the same dense grid as in 
the first phase to acquire the UWB waveform in 
^aaj = 105 distinct vertices. The positions of the 
transmit antenna in the two phases were (0.26, 2.97) 
m and (0.26, 3.97) m, respectively, if the coordi- 
nate system shown in Fig. [T] is employed. The 
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Figure 1 . Map of the enviroment of our experimental campaign. 



distance between the two antennas varied between 
1 m and about 5 m; larger distances were not 
taken into consideration since we were interested in 
indoor ranging only. It is important that to note that 
the choice of the measurement scenarios described 
above is motivated by the fact that UWB signals 
experienced similar propagation in both phases. 

For each position of the receive antenna up to 
25 realizations of the received UWB signal were 
acquired (each one lasting Tacq = 110 ns), in order 
to filter out the effects of the measurement noise 
and the TOA estimate outliers in the postprocessing 
phase. The analysis of the multiple waveforms re- 
ferring to each link evidenced the presence of few 
time-variant artifacts (mostly due to the movements 
of the people behind the radios); as a matter of 
fact, the channel in our measurement campaign 
can be deemed quasi-static and is characterized by 
reflections mainly due to the walls, the ceiling and 
the floor of the room (or rooms) hosting the radio 
devices. 

Finally, for each acquired waveform, a TOA 
estimate (evaluated by the PulsON220 devices) and 
the actual transmitter-receiver distance (evaluated by 
means of a metric tape with an accuracy better than 
1 cm) were also stored in the database. The TOA 
estimates have been used to provide a common time 
frame to all the acquired waveforms; this has made 
possible the estimation of the mean excess delay 
and of other signal statistics in the signal processing 
phase. 
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B. Statistical modeling of signal features 

In each of the two measurement phases de- 
scribed in the previous Paragraph, the procedure 
for acquiring UWB waveforms worked as follows. 
First, the UWB radios carried out their handshaking 
procedure for achieving mutual time synchroniza- 
tion; then, the receiver sampled the incoming wave- 
forms and generated a TOA estimate for each of 
them; let r,^°^(t) (rf^°^(t)) and r[°^ = TOA{rf'{t)) 
(^NLos A TOA(rf°^(t))) denote the received wave- 
form and the corresponding TOA estimate (here the 
function TOA(-) represents the internal procedure 
adopted by the PulsON220 devices, based on energy 
thresholds and a go-back technique, to estimate the 
TOA; see [9] for further details) acquired over the 
i-th link in LOS (NLOS) conditions. The modeQ 
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was adopted for the measured TOA (for both LOS 
and NLOS conditions). Here, cq is the speed of 
light, di denotes the distance between the transmitter 
and the receiver, 6j is the NLOS bias (in seconds) 
affecting the TOA measurement (the values taken 
on by this random parameter are always positive 
for NLOS links and null for LOS hnkg and 
Wi ~ A/'(0,cr^j) is the measurement noise; in 
addition, the expression a^- = 70"^ (if is adopted 
for the variance of the measurement noise, where 7 
is a parameter depending on both the specific TOA 
estimator employed in the ranging measurements 
and on various parameters of the physical layer, 
and (3 is the path-loss exponent (a known and fixed 
value is assumed for this parameter in both LOS 
and NLOS conditions [4J). 

In the following we focus on the problem of 
estimating the bias hi (affecting the TOA estimate 
Ti ([T])) from a set of A^/ = 6 different "features" 
{xij, j = 0, 1, ..., 5} extracted from the received 
waveform ri(t). In particular, like in []2l, the fol- 
lowing features were evaluated for the set of the 
received waveforms: 



3) the delay spread Xi 2 

rm,?^^dt; 

4) the energy Xj,3 = Ei = J^ \ri{t)\^dt; 

5) the rise time Xj4 = trise,i 
min{t : |rj(t)|/maX( |rj(t)| > 0.9} — 
{t: \ri{t)\/ma.xt\ri{t)\ > 0.1}; 

6) the kurtosis Xjs = Kj = 
^/r (kiWI -Ai|r|) dt, where fi\r\ = ^ 

\r\ 

Jr^\ri{t)\dt, af^i = ^ Jrpi\ri{t)\ - iJi\r\fdt and 

T denotes the observation time. 
In the following the set {rj(t)} of received wave- 
forms is modelled as a random process, so that the 
above mentioned features form a set of correlated 
random variables; in addition, all of them are sta- 
tistically correlated with the TOA bias. The last 
consideration is confirmed by the numerical results 
of Table [1} which lists the absolute values of the 
correlation coefficients of the previously described 
features with the estimated TOA bias for both the 
LOS and the NLOS scenarios. From these results it 
can be easily inferred that not all the considered 
features are equally useful to estimate the TOA 
bias. For this reason and to simplify our statisti- 
cal analysis, we restricted the group of considered 
features to the set [xi^, Xj_i, Xi^2}, which collects 
the parameters exhibiting a strong correlation with 
the bias in the NLOS scenario. 

Note that, in principle, the bias hi is not influ- 
enced by the distance di, since it depends only on 
the thickness of the walls (or of other obstacles) 
encountered by the transmitted signal during its 
propagation; this is not true, however, for the above 
mentioned triple of signal features (see [1] for 
further details). Generally speaking, it is useful to 
derive a TOA bias estimator which is not influenced 
by the transmitter-receiver distance di. In the at- 
tempt of removing the dependence of the features 
{xifi, Xi^i, Xi^2} on the link distance, we developed 
the models 



1) the maximum signal amplitude Xj.o 
maxi \ri{t)\; 
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'The superscripts LOS and NLOS are not explicitly indicated in the 
following expressions, when not strictly needed, to ease the reading. 

^In this case the probability density function (pdf) of bi is /i,(6) = 
S{b). 
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Correlation 
with bi 


Xi,0 


Xi^l 


Xi,2 


Xi,3 


Xi,4 


Xi,5 


d. 


NLOS case 


0.795 


0.852 


0.894 


0.641 


0.454 


0.644 


0.609 


LOS case 


0.602 


0.586 


0.129 


0.666 


0.586 


0.119 


0.629 



Table I 
Absolute VALUES of correlation coefficients between bi ANDxij, with j = 0, ....,Nf 



on the basis of our experimental results (and, in 
particular, on the basis of the estimates of the 
joint pdf's {/;,,^.„, fh,xi, fb,:^^} referring to the three 
possible couples {bi, Xij) with j = 0, 1, 2); here 
r^ax' '^m ^^'^ '^Ts ^6 random variable^ whereas r^^ 
and r™„^ are deterministic parameters having known 
values. Given these models, the vector of distance- 
independent features x = [rma^.,r™, r^] can be 
evaluated from its distance-dependent counterpart 
X = [rmax-, Tm, ^ds] and can be used in place of it. 
A complete statistical characterization of the es- 
timated bias and of the 3 related signal features we 
consider is provided by the joint pdf fb,rma.,rm,rA-) 
or, equivalently, by the joint pdf /fe,rO,,,,r-,r™(-); 
these functions have been estimated from the ac- 
quired data, but cannot be plotted because of the 
large dimensionality of their domains; for this rea- 
son, we analysed some pdf's deriving from their 
marginalization. For instance. Fig. |2] shows the 
estimated joint pdf's for the estimated bias bi and 
measured delay spread for the NLOS and the LOS 
scenarios; these results deserve the following com- 
ments: 

1) a significant (limited) correlation between 
these parameters is found in the NLOS (LOS) 
case; 

2) the null region exhibited by the estimated 
pdf's is due to the fact that the TOA bias can- 
not take on values in the interval [0, t^aii/co], 
where t^aii is the thickness of the wall ob- 
structing the direct path; 

3) large values of the TOA bias are unlikely 
since they are associated with small incidence 
angles of the transmitted signal on the ob- 
structing wall. 

Finally, it is important to point out that our exper- 
imental data were also processed to evaluate the 
correlation coefficient between the estimated bias 
bi and each of the distance-independent parameters 

^In the rest of the document, the subscript i has been omitted for 
simplicity when not strictly necessary. 
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Figure 2. Estimated joint pdf of the estimated TOA bias and 
the delay spread in NLOS conditions (above) and LOS conditions 
(below). Note that the NLOS PDF has been scaled by a factor 3 to 
improve its visualization. 
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Figure 3. Estimated joint pdf of the estimated TOA bias and 
the r^£,2, random variable in NLOS conditions (above) and LOS 
conditions (below). Note that the NLOS PDF has been scaled by a 
factor 2 to improve its visualization. 



r^ax' '^m ^iid r^ (see Table |ll] and Fig. [3} which 
represents the joint pdf of the estimated bias 6j and 
the parameter r^„^); our results evidence that the 
last parameters are less correlated with bi than their 
distance dependent counterparts. For this reason, we 
decided to take into consideration also the distance- 
dependent features {xj 0, Xi^i, Xi^2] for localization 
purposes. 
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Correlation 
with bi 


Xifi 

' max,t 






NLOS case 


0.321 


0.649 


0.561 


LOS case 


0.483 


0.370 


0.490 



Table II 
Absolute values of correlation coefficients between 

THE ESTIMATED BIAS it AND EACH OF THE PARAMETERS rj^^^,, 

T^ AND r^. Both LOS and NLOS conditions are 

CONSIDERED. 



III. Localization algorithms 
A. Introduction 

In this Section we develop various algorithms 
for two-dimensional localization in a UWB network 
composed by A^^ anchors with known positions 
z^ 4 [x^^ y^f e M^ i = 0,...,Na - 1 and by a 
single node (dubbed mobile station, MS, in the fol- 



lowing) with unknown position 6 = [x, y] G M . 
Any couple of the given (Na+l) devices can operate 
in a LOS (NLOS) condition with probability Plos 
(1 — Plos)- The localization algorithms described 
below try to mitigate the effects of the NLOS 
bias error and aim at generating an estimate 
of 6 minimizing the mean square error (MSB) 



6 — 6\\ k It is also important to point out 
that localization algorithms developed for NLOS 
scenarios usually consist of two steps. In fact, first 
the NLOS bias is estimated for each involved link 
and is used to remove the bias contribution in the 
acquired data; then the new data set is processed 
by a least-square (LS) procedure generating an 
estimate of (e.g., see @, [H, S, [H). A different 
approach, involving implicit estimation of the bias 
for each link, is adopted in the following; this 
approach is motivated by the fact that the estimation 
of the bias for the i-th link can benefit from the 
information acquired from the other (A^^ — 1) links; 
note that in [2J and [|61 bias mitigation performed in 
a link-by-link fashion is exploited to assign a weight 
in a weighted least-square (WLS) step, but such an 
approach is heuristic. 



B. Maximum Likelihood Estimation 

If the links between the MS and the Na different 
anchors are assumed mutually independent, the ML 
estimation strategy of the unknown MS position 0, 
given a TOA estimate and a set of additional signal 



features for each link, can be formulated as 

Na-l 

^ = argmaxln TT /^,x(i^.,Xi;^). (5) 

^-^ 



i=0 



Here, 



^ n~AT 



\x 



yy 



denotes the MS trial 



position, Ti and Xj are the TOA and A^/-dimensional 
signal vector collecting the received signal features 
acquired for the i-th link and /t-x(-,-;^) is the 
joint pdf of the TOA and the vector of features 
parameterized by the trial position 6. As already 



discussed in Paragraph |II-B[ various options for x 
(and, consequently, for fr,:i<.{-, ■; 0)) are possible; in 
the following Paragraphs the impact of such options 
on the ML strategy (|5]) is investigated. 

1) ML estimation strategy for different sets of ob- 
served data: In this Paragraph two different options 
are considered for the set of features processed by 
the ML strategy. 

a) Option A: In this case it is assumed that 
the vector of features referring to the i-th link is 
Xi = [rmaxA, TmA, Tds,i\'^; tMs choicc is motivated by 
the large correlation between these random variables 



and the link bias b (see Paragraph II-B). Then, 



the joint PDF fr^yi{Ti,yii]6) appearing in the ML 
strategy (|5]) can be expressed as (see Eq. ([T])): 



/r,x(T-i, Xi; 0) = (/b,x ® /^) I Ti - 

where (see Eq. Q, ([3]), (Q) 



''"maxA ~r T 



MO) 



T„ 



dm 

Co 



'^ds,i T^s 



(6) 



d^io) di{e) 



(7) 

di{6) is the distance between the i-th anchor and the 
MS trial position and /^ x /«, denotes the convolu- 
tion between the joint pdf /;, x(-) and the observation 
noise pdf f^iw) = (27ra2)-V2exp(-^V(2a2)). 

Note that the shape of the function fb,x{-) under 
the hypothesis of LOS conditions (i^Los event) is 
substantially different from that found in NLOS 
conditions (i^NLos event) [^; in both cases this func- 
tion was estimated applying the procedure described 
in [m to the data collected in our measurement 
campaign. This led to two distinct multidimensional 
histograms, which approximate the required pdf's 
fb,iib,x,\HLos) and /^.x (^XjI-^nlos) with a certain 
accuracy depending on: a) the amount of acquired 
data; b) the sizes Ah, Armax, ^^m and Ar^^ of 
the quantization bins adopted in the generation of 
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the histograms. Note that these sizes need to be 
accurately selected, since large bins imply a coarse 
approximation of pdf s, whereas excessively small 
bins require a huge amount of data. 

Given an estimate of the above mentioned couple 
of pdf 's, the pdf /6,x(-) can be evaluated as 

/fe,x {b, Xi) = Pms/b.x {b, Xi|ifLOs) + 

(l-ilos)/fe,x(&,i^|i/NLOs) (8) 

if the probabilities Plos = Prji^Los} and Pnlos = 
Prji^NLos} = 1 — Prj-f/^Los} are available, or as 

fb,k (b, Xi) = -fb,ic {b, Xili^Los) + ^fb,k (b, Xili^NLos) 

(9) 
if no a priori information about the LOS/NLOS 
conditions are available. 

In principle, evaluating the joint PDF 
fT,yLiTi,^i'-iG) of the ML strategy entails the 
computation of the convolution appearing in (|6]); 
this has not been done in our work, since it is 
numerically complicated and a form of implicit 
smoothing is already included in the generation of 
the above mentioned histograms extracted from the 
available measurements. 

b) Option B: In this case the set of features 
employed in ML estimation consists of a single 
element, namely the delay spread (which exhibits 
the largest correlation with the NLOS bias; see 
Table |lf, so that Xj = Tds,i and the pdf /r,x(rj, x^; 0) 
of ([5]) becomes (see ([6j)) 



frA'^i^^uO) = (/fe,T-®/«;) 



dm 

Co 



Tds 



'ds 



dAO) 



(10) 

Like in the previous case, the pdf fb,T^{-) was 
estimated in the LOS and NLOS scenarios (see 
Fig. |2]) from the data acquired in our measurement 
campaign. Note that this option leads to a ML 
localization algorithm which is substantially simpler 
than that proposed in the analysis of option A. 

2) Parameterization of the observations: The 
ML strategies developed above are based on the 
joint pdf's hA^^,ri^,rz{-) and /6,r™(-). which refer 
to a set of distance-independent parameters. Since 
the parameters r^^^, r™ and r^ exhibit a lower 
correlation with b than their distance-dependent 
counterparts Vmax, Tm, T~ds (see Tables |I] and |II]), the 
use of the joint PDF 



fT,^{ri,yii]e) = (/6,x®/« 



r,-*^.x,|(n) 



in place of ([6]) has been also investigated (similar 
comments apply to ([10])). 

3) Estimation of joint pdf's: As already ex- 
plained above, the joint pdf's involved in the pro- 
posed ML localization strategies can be easily es- 
timated from the acquired data using a simple 
procedure based on dividing the space of observed 
data in a set of bins of proper size. Such a pro- 
cedure generates an histogram, which, unluckily, 
entails poor localization performance if employed 
as it is, because of the relatively small number 
of bins (adopted to avoid empty bins). To miti- 
gate this problem, interpolation followed by low- 
pass filtering can be applied to raw experimental 
histograms. In particular, we found out that cubic 
spline interpolation followed by moving average 
filtering (with a window size equal to about 1/15 
of the size of the histogram) provides good results 
(an example of the resulting pdf is shown in Fig. 
|2]); the localization strategies adopting this approach 
are dubbed interpolated-histogram estimators in the 
following. 

An alternative to this approach is represented by 
fitting the histograms with analytical functions (e.g., 
polynomials). On the one hand, this solution offers 
some advantages with respect to the previous one, 
since it requires less memory (only the values of the 
coefficients of the fitting equations need to be saved) 
and provides better accuracy (the resulting pdf's are 
smoother). On the other hand, the use of (multi- 
dimensional) fitting raises a number of problems. In 
fact, global performance figures, like the root mean 
square error^ (RMSE) are certainly useful, but do 
not account for the fact that fitting errors where a 
pdf takes on small values are less critical than those 
in points where the same pdf takes on large values. 
In addition, if the dimensionality of the pdf domain 
is larger than two, an analytical function generated 
by a fitting procedure cannot be easily represented. 

In our statistical analysis, the use of polyno- 
mial fitting was investigated. A polynomial degree 
Npoi = 8 was selected, since this choice ensures that 
the resulting RMSE is smaller than t fitting = S-IO^"^ 
in both LOS and NLOS scenarios. In the follow- 
ing localization algorithms employing pdf's gener- 
ated by polynomial fitting will be dubbed fitted- 
histogram estimators. 

'^This parameter is evaluated as the root mean square of the differ- 
ence between the values of the original experimental histogram and 
the corresponding values generated by a fitting analytical function. 
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Assume 

H = H,„, A = 



Form the 
decision 
variables: 




Compute the posterior probabihties: 

P,=\f,Jb,^\HLos)dh 

P,., = lA.,(b,Xi I H,„,)dx, j = 0,1,...,JV, - 1 

9,„=lA,(b,'',\H,,os)dx, J = 0,l,...,N, -1 



H = H„ 



H = H„ 



Recompute the bias 
estimate: 

b - argmax/^^(A,x | H - H) 



Figure 4. Flow diagram of the proposed iterative estimator. 



C. Iterative Estimation 

Recently, an iterative estimator of both the chan- 
nel state (i.e., LOS or NLOS conditions) and the 
NLOS bias h has been proposed in [4J. Such an 
estimator relies on the availability of univariate 
statistical models for the estimates of TOA, RSS and 
RDS; it is shown, however, that estimation accuracy 
is improved if the delay-spread is exploited as the 
only discriminant feature. 

In the following, a modified version of this itera- 
tive algorithm (dubbed iterative estimator), employ- 
ing the joint pdf's described above and extracted 
from our experimental database, is proposed. This 
algorithm is summarized by the flow diagram shown 
in Fig. |4] and operates as follows. In its first iteration 
it starts with a null estimate of the NLOS bias 
and, for each link, generates a set of posterior 
probabilities integrating the pdf's estimated in LOS 
and NLOS conditions. Then, such probabilities are 
employed to compute a couple of decision vari- 
ables on the basis of which a decision H on the 
LOS/NLOS conditions is taken. Finally, such a 
decision is exploited to compute a new bias estimate 
which is used to start the next iteration. Note that 
this algorithm operates in a link-by-link fashion and 
that in the diagram of Fig. |4] the dependence on the 
link index has been omitted to ease the reading. 

IV. Numerical results 

Our experimental database has also been ex- 
ploited to assess the RMSE performance of the 



proposed algorithms for localization and NLOS 
bias mitigation schemes via computer simula- 
tions. In our simulations the MS coordinates are 
always (0;0); then, following [2J, for the i- 
th anchor a received waveform from either the 
LOS database (with probability Plos) or from the 
NLOS database (with probability Pnlos) was drawn 
randomly and was associated with the position 
z^ = ((iisin(27r^); di cos(27r^) j, where di is the 
distance measured for the selected waveform. Note 
that these waveforms already include the experimen- 
tal noise and thus no simulated noise was imposed 
on the waveforms (so that the signal-to-noise ratio 
is the experimental one). Finally, the parameter Na 
has been set to 3 (worst case which still theoretically 
allows unambiguous localization). 

The RMSE performance of the following algo- 
rithms has been evaluated: 



1) 



2) 



3) 



4) 



5) 



LS - A standard LS estimator for LOS en- 
vironments; the estimation strategy can be 
expressed as = argming^.^Q" (coTj — 

di{e)f. 

VE - A LS estimator exploiting TOA mea- 
surements corrected by the algorithm pro- 
posed by Venkatesh and Buehrer in [4J; this 
algorithm relies on a statistical modeling of 
the propagation environment based on our 
experimental database. 

ML-4D - A ML estimator based on (|9j) and 
employing an interpolated histogram for the 
evaluation of its likelihood function referring 
to a distance-dependent parameterization. 
ML-2D - A ML estimator based on (fTOb and 



employing an interpolated histogram for the 
evaluation of its likelihood function referring 
to a distance-dependent parameterization. 
ML-2D-ID - A ML estimator based on (fTOb 



6) 



7) 



and employing an interpolated histogram for 
the evaluation of its likelihood function refer- 
ring to a distance-independent parameteriza- 
tion. 

ML-4D-F - A ML estimator based on ^ and 
employing a fitted histogram for the evalua- 
tion of its likelihood function referring to a 
distance -dependent parameterization. 
ML-4D-IT - A LS estimator based on (|9j) 
and exploiting TOA measurements corrected 
by the modified iterative algorithm illustrated 
in Section UlLCl 
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8) ML-2D-IT - A LS estimator based on (10) 



and exploiting TOA measurements corrected 
by the modified iterative algorithm illustrated 
in Section UlLCl 

In estimating the RMSE performance of the ML 
algorithms listed above the likelihood functions 
were always evaluated at the vertices of a square 
grid characterized by a step size equal to 10 mm. 
Some numerical results are compared in Fig. |5} 
which illustrates the RMSE performance versus the 
probability Plos- These results evidence that: 

1) Different estimators can provide significantly 
different accuracies for distinct values of Plos, 
when this probability is not close to unity. 

2) The simple LS algorithm is outperformed by 
all the other algorithms when Plos < 0.9; this 
is due to the fact that this strategy does not 
try to mitigate the NLOS bias. 

3) The VE algorithm performs well at the cost 
of a reasonable complexity, but offers limited 
bias mitigation when Plos = 0; in this case the 
ML-2D, ML-2D-ID and ML-4D-F estimators 
perform much better. 

4) The exploitation of a large set of received 
signal features does not necessarily allow to 
achieve better accuracy than a subset of them 
(see the curves referring to ML-2D and ML- 
4D estimators); this is due to the fact that the 
correlation between the different couples of 
extracted features is typically large, so that 
they provide strongly correlated information 
about the NLOS bias. 

5) Distance-dependent parameterization provides 
better accuracy (see the curves referring to the 
ML-2D and ML-2D-ID estimators); this can 
be related to the fact that r^ is less corre- 
lated with the NLOS bias than its distance- 
dependent counterpart t^s- 

6) The ML-4D-F estimator performs better than 
the ML-4D estimator in NLOS conditions; 
this means that the use of fitted histograms 
entails an improvement in localization accu- 
racy. 

Finally, Fig. [6] shows the cumulative density func- 
tion (CDF) of the RMSE localization error charac- 
terizing the LS, VE, ML-2D and ML-2D-IT algo- 
rithms when Plos = 0.2. These results evidence that: 
a) the RMSE localization error of the considered 
algorithms remains below 0.5 m in average the 
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Figure 5. RMSE performance versus Plos offered by various 
localization aleorithms. 
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Figure 6. CDFs of the RMSE characterizing various localization 
algorithms when Plos = 0.2. 



80% of the cases; b) the ML-2D technique is more 
accurate than the other algorithms, in terms of 
RMSE, for small NLOS errors. 



V. Conclusions 

In this paper various UWB localization tech- 
niques processing multiple features extracted from 
the received signal to mitigate the problem of NLOS 
bias have been described and their accuracy has 
been assessed exploiting the experimental data ac- 
quired in an measurement campaign. Our results 
evidence that: a) a restricted set of features has to be 
employed; b) the use of distance-dependent features 
and of fitted histograms provide better performance 
than that offered by distance-independent features 
and interpolated histograms. 
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